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We re-examine and correct an earlier derivation of the dis- 
tribution of the Wigner phase delay time for wave reflection 
from a long one-dimensional disordered conductor treated in 
the continuum limit. We then numerically compare the dis- 
tributions of the Wigner phase delay time and the dwell time, 
the latter being obtained by the use of an infinitesimal imagi- 
nary potential as a clock, and investigate the effects of strong 
disorder and a periodic (discrete) lattice background. We find 
that the two distributions coincide even for strong disorder, 
but only for energies well away from the band-edges. 



I. INTRODUCTION 

The delay time associated with potential scattering is 
one of the important quantities related efficiently to the 
dynamical aspect of scattering in quantum mechanics. 
One of the common measures for this quantity is the 
Wigner phase {(j)) delay time — fi(dcj)/dE)\ [|], which 
essentially entails following a fiducial feature such as the 
peak of the wavepacket as it traverses the scattering re- 
gion. This procedure is, however, rendered meaningless 
under conditions of strong distortion of the wavepacket 
by the scattering potential Further, there is the 

problem of indentifying the position of the particle with 
the peak of the wavepacket. Several researchers have 
made other proposals for identifying a physically mean- 
ingful timescale of interaction of the particle with the 
scattering potential (See for recent reviews These 
include the quantum clocks that utilize the co-evolution, 
in a locally applied infinitesimal field / potential, of 
an extra degree of freedom (such as the spin |^) at- 
tached to the traversing particle. Even these propos- 
als are not completely free from problems |§j|J^,^. One 
related quantity that has, however, remained uncontro- 
versial is the dwell time obtained with the 'non-Unitary' 
clock [||,^[j9j involving absorption/amplification due to a 
locally applied infinitesimal imaginary potential [Vi) 
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for the case of total reflection (|i?p = 1). fn this case, 
it also turns out to be the average Smith dwell time 
Td = i^/j) / l^pc'a^, where j is the incoming flux in the 
steady state situation, and ip is the wavefunction in the 
scattering region of interest. In this paper, we will con- 
sider the Wigner phase delay time and the dwell time, 
given by the non-Unitary clock, for total reflection from 
a long one-dimensional disordered medium. These times 



are, however, not self-averaging and one must have their 
full probability distribution over a statistical ensemble of 
random samples. 

The distribution of these times for the random media 
has been investigated recently by several workers [p[-p2[ . 
A delay time distribution that appears universal for wave 
reflection from a long one-dimensional (one-channel) ran- 
dom system was derived recently by Texier and Comtet 
]l6t in the limit of high energy(i<^) and weak disorder as 



(2) 



where a — 4(A^fc)^^, is the strength of the dis- 
order (see Eq. (6)), and the dimensionless delay time 
T = ET^/h. This was later confirmed by Ossipov et 
al. [p^ for a discrete random chain. 

Earlier, we had derived the distribution of the dwell 
time for total reflection, i.e., in the insulating limit, using 
the non-Unitary clock for both passive as well as active 
(absorbing or amplifying) one-dimensional random con- 
tinuous media The dwell time distribution obtained 
by us, under the condition of a random phase approxi- 
mation which is valid for high energy and weak scatter- 
ing, coincided exactly with the delay time distribution of 
Texier and Comtet. However, an earlier calculation by 
Jayannavar et al. pO| ] for the Wigner phase delay time 
had obtained a slightly different form from the distribu- 
tion and it has been speculated that this discrepancy 
may well be due to the continuum model used by them. 
In this paper, we first re-examine and correct the calcu- 
lation of Jayannavar et al. for the Wigner phase delay 
time distribution. We find that the discrepancy noted 
above arises actually from an inconsistency of the approx- 
imations made within the random phase approximation 
(RPA), and that when the approximations are carried out 
consistently, their expression reduces to the universal dis- 
tribution as in Eq. (^ . We then examine numerically the 
distribution of the delay times and the dwell times for 
strong disorder using the transfer matrix method for a 
one-dimensional disordered chain with a one-band tight 
binding Hamiltonian. We find that the distribution of the 
Wigner phase delay time and the dwell time, clocked by 
the non-Unitary clock, agree exactly even in the strong 
disorder regime, but for energies far away from the band- 
edge. We further examine the effect of a periodic lattice 
on the delay time by varying the energy within the band. 
We find that for strong disorder, and for energies close 
to a band-edge, the Wigner phase delay time distribution 
differs considerably from that of the dwell time given by 
the Non-Unitary clock. The Wigner phase delay time 
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can even become negative under such conditions. The 
dwell time, however, remains positive as it must for total 
reflection (|i?| = 1). 



where the dimensionless length I = L/lc = 1/2 A^fc^L. 
In the limit of large lengths, Z 3> 1, the distribution sat- 
urates and we can set dP/dl = 0. Hence, we obtain the 
solution llOl 



II. DISTRIBUTION OF THE WIGNER PHASE 
DELAY TIME FOR TOTAL REFLECTION IN 
THE RPA 

Here we will re-examine the earlier calculation of 
Jayannavar et al. (l^Jl^ for the distribution of the 
Wigner phase delay time for total reflection from a one- 
dimensional (one-channel) disordered continuum. Again, 
we will begin with the invariant imbedding equation for 
the reflection amplitude R{L) = \/r{L) exp[i0(L)] given 
by (in the notations of Ref. pX(|) 



dR{L) 
dL 



2ikR{L) + j7^AL)[l + RiL)f, (3) 



where rir{L) = —Vr{L)/E is the normalized fluctuating 
potential. In the limit of large lengths (L ^ Ic, the lo- 
calization length), the reflection becomes approximately 
total {r{L) ~ 1), and Eq.(^ yields an equation for the 
reflection phase (t>{L) as 



del) 
dL 
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The equation for the phase delay time = h{d(j)/dE) = 
l/cg(d<p/dk) (where Cg is the group velocity), is obtained 
by differentiating the above equation for <j) with respect 
to k: 

^ = — [2 + 77^(i)(l + cos0-/cCgr0sin0)]. (5) 

Cl-LJ Cg 

As before, we will assume the random refractive index 
rir{L) to be a Gaussian white noise with a zero mean, 
i.e., 



Poo(ti) = 



(eW2_i)(i + r2) 



(8) 



where A — 8/\/3A^fc and the dimensionless time ti = 
CgkT^/V3. This expression does yield the r-^ ^ universal 
tail behaviour for ri — > cx), but differs from the distribu- 
tion of dwell times given by Eq.(|^) at short times. The 
main difference appears at t = 0, where this expression 
for Poo (t) yields a flnite value in contrast to Eq. ^ which 
gives P^{t) = because of the essential singularity at 

T=0. 

This difference is readily traced to the fact that the 
RPA is good only in the high energy, weak disorder 
limit. Indeed, if we consistently take the high energy, 
weak disorder Hmit in Eq. (0) [or in Eq. (||)], i.e., by de- 
manding Cg/c^ oo and jcg — > with the product 
(A^/cg)(CgA:^) = 4/q; a constant, we obtain the solution 
as 



Poo(r) 



a 



exp( ), 



(9) 



where t = ET^jti. This is exactly the full imiversal dis- 
tribution of delay times obtained in Eq. (||) for the 
case of a free electron with Cg = hk/m, thus reconfirming 
again the universal delay time distribution. We note that 
the above approximations have to be carried out consis- 
tently specially for a large group velocity Cg , which is the 
case for energies far away from the band edges. This also 
suggests that the condition of weak disorder A^/c ^ 1 for 
the one-parameter scaling, which assumes a uniform dis- 
tribution of the phase (RPA), may have to be modified 
to t^kjicgjc^) <C 1, where is the phase velocity. 



(7],(L)) = 0, ML)7jr{L')) = A^6iL - L'). 



(6) 



Using the Novikov theorem p^ , we can now set up a 
Fokker-Planck equation for the joint probability distri- 
bution function P{T^,(j)\L) over the ensemble of rir{L). 
However, we will be interested in the marginal proba- 
bility distribution P{T^]L), which can be obtained by 
integrating over the phase angle 4>. To this end, we 
make the random phase approximation (RPA) and set 
P{Tcj,, (p; L) = P{T^; L) /2Tr, i.e., assume a factored out 
uniform distribution over the phase angle 4>. The RPA is 
a good approximation for high energy and weak disorder 
p4| . We obtain the equation for P{T^] L) as 
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III. STRONG DISORDER AND A PERIODIC 
BACKGROUND: NUMERICAL RESULTS 

The probability distribution of dwell times in Ref. ||] 
was derived for a continuum model in the limit of weak 
disorder and high energy when the RPA is valid. In this 
Section, we will examine these limitations numerically. 
In particular, we investigate the distributions of dwell 
times for the case of strong disorder and compare the 
distributions of the dwell times and the Wigner phase 
delay times. We will simulate a disordered lattice, in- 
stead of a continuum. The underlying lattice will also 
provide a discrete periodic background in the system, as 
distinct from a uniform continuum, whose effect on the 
delay times will be investigated. 

In order to go beyond the RPA, we will use the transfer 
matrix method involving the products of random trans- 
fer matrices 125(1 to simulate the one-dimensional random 
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medium using the one-band tight binding Hamiltonian 
(the Anderson Hamiltonian) with diagonal disorder . 
The Hamiltonian describing the motion of a particle on 
the random lattice can be written as 

H = ^[e„|n) (n| + V{\n) + + + {n\)] (10) 

n 

where |n), e„ and V denote, respectively, the non- 
degenerate Wannier orbital at the v}^ site, the site energy 
at the n*^ site and the hopping matrix element connect- 
ing the nearest neighbours separated by a unit lattice 
spacing. The site energies e„ can be written explicitly in 
the form of — iry, with the real parts of the site energies 
assumed to be independent random variables distributed 
uniformly over the range [—W/2^W/2] for 1 < n < iV 
and zero otherwise. This is so that the disordered chain 
of N sites is imbedded in an otherwise infinitely long 
ordered lattice. The imaginary part in the site energy 
(—17/) makes the Hamiltonian non-Hermitian and causes 
the particles to be formally coherently absorbed or emit- 
ted depending on the sign of 7], which is taken to be con- 
stant and non-zero (though infinitesimally small) over the 
disordered segment 1 < rt < iV and zero elsewhere. Since 
all the energies can be scaled with respect to V , we will 
set V to unity. 

The reflection [R) and the transmission(r) amplitudes 
can now be calculated using the transfer matrix method 
p5| . In order to calculate the Wigner phase delay time, 
the reflection amplitude R{E) = \/r{E) exp[— i0(£')] is 
computed at two slightly differing values of the incident 
wave energy, E = Eq and E — Eq + 6E, for a conser- 
vative chain (r/ = 0). The Wigner phase delay time is 
then calculated as = h{d4>/dE) = h[(j){Eo + 5E) - 
(^{Eq)]/ 5E. Similarly, to calculate the dwell time by ap- 
plying the imaginary potential, the reflection amplitude 
is computed at two values of the imaginary site energy 
(?7 = and rj = Srj). Now the dwell time is given by 
Td = h/2{d\R\ydrj) = h/2[\R{E,7j = St^)]^ - \R{E,r, = 
0)p]/57y. Typically the values of SE and 6r] are 10~^ and 
the stability of the results have been checked for their 
choice within the range 10~^ < SE, St] < 10"''. (We will 
deal with the delay dwell time in a dimensionless form 
by multiplying it by V and setting fi ~ 1.) For the cal- 
culation of the averages and the distributions, we have 
typically used 10^ configurations of the disorder. We will 
present our results for a long sample {L :s> Ic, i.e., lengths 
much greater than a localization length). 

We will first examine the case of wave energies far away 
from a band edge {E = 0.0, 1.0). In Fig. ??, we show 
the distribution of the Wigner phase delay time Tw and 
the dwell time for reflection from a long sample for dif- 
ferent values of the disorder strengths {W = 0.1, 2). For 
weak disorder {W — 0.1), the distributions are identical 
to each other and also correspond exactly to the univer- 
sal distribution given by Eq.(^. It should be noted that 
the RPA is not valid for exactly the band centre {E = 0) 
due to a well-known anomoly, although it is valid for a 



generic value of energy within the band |2^]. The two 
distributions also coincide for higher disorder strengths. 
psf It is interesting to note that Eq. (||) still describes the 
distribution reasonably well for moderately large disorder 
{W = 2.0, see Fig. ??b), though the RPA under which 
the expression was derived is not valid for these cases. 
The case of E = 1 shows similar behaviour, though the 
peak occurs at a different value, reflecting the smaller 
group velocity. In Fig. ??d, we plot the distributions of 
dwell and delay times for cases of symmetric disorder (the 
random site energy is chosen from the range symmet- 
ric about zero [—W/2,W/2]) and asymmetric disorder 
(the site energy can only be positive [0, W] or negative 
[—W, 0]). The distributions for the positive and the neg- 
ative one-sided, asymmetric disorder appear to be the 
same, regardless of the sign as expected, of course (we 
have included these two only as a check on our numer- 
ics). These, however, are different from the distribution 
for the symmetric case. The contribution of the prompt 
part of the reflection arising from the average potential 
mismatch at the boundary is clearly seen for the cases of 
asymmetric disorder in that the peak of the distribution 
occurs slightly earlier, and there is more weight at early 
times. The asymmetric disorder model changes on an 
average the lattice potential locally only over the disor- 
dered segment (sample) and not globally over the sample 
and the leads. This gives rise to a potential mismatch at 
the ends of the disordered chain and cannot be absorbed 
by a mere shift of the incident energy. More formally, 
in the case of the symmetric disorder the average of the 
S-matrix, (S) — 0, while in the case of the asymmetric 
disorder, (5) ^ 0. 

Now, we will examine the case of wave energies close 
to the edge of the band {E = 1.9, 1.99). In Fig. (??), we 
show distributions of the delay time and dwell times. For 
the case of weak disorder, again the Wigner delay time 
distributions and the dwell time distribution coincide. 
There is, however, considerable discrepancy from Eq.(^, 
as can be seen. We have explicitly verifled that the RPA 
is valid for this case of weak disorder by calculating the 
distribution of the phase. Thus, the discrepancy cannot 
be an artefact of the RPA. The most probable reason, 
perhaps, is that near the band edge, the wave does not 
penetrate deep enough to fully sample the randomness, 
before getting reflected promptly. This would call into 
question the factorization of the joint probability distri- 
bution of the phase and its energy derivative, particularly 
for short times. 

For intermediate and strong disorder {W — 1, 2, 3), a 
more interesting effect occurs. The two distributions, i.e., 
the Wigner delay time distribution and the dwell time 
distributions no longer coincide. The difference between 
the distributions increase with the disorder strength and 
with their proximity to the band-edge. The Wigner delay 
time distribution appears quite different from the univer- 
sal distribution at i? = 0. Near the band-edge, in fact, for 
E = 1.99 and W = 1, the Wigner delay time distribution 
is non-zero for even negative times. This is, presumably. 
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due to the strong deformation of the wavepacket caused 
by the strong dispersion near the band-edge. The dweU 
time distribution given by the 'non-Unitary' clock, how- 
ever, remains non-zero only for positive times. We also 
note that the Universal tail at long times (r — > oo) 
remains unaffected. 



IV. CONCLUSIONS 

In conclusion, we have studied the distribution of the 
delay and the dwell times for reflection from a disordered 
medium in the limit of total reflection {\R\ = 1). We have 
revisited the original calculation of Jayannavar et al. | [Tc| ] 
for the distribution of the Wigner delay time. We show 
that, by taking the high-energy limit consistently within 
the RPA, the correct universal distribution of delay times 
is reproduced. In the course of the derivation, we note 
that the single-parameter scaling ansatz for the RPA 
seems consistent under the condition A'^k/{cg/c^) <^ 1 
(A^ -the disorder strength, Cg-the group velocity and 
C0-the phase velocity), instead of A^fc ^ 1 which does 
not account for the effects of the group velocit^^ This 
is in accord with the recent results of Ref. |l^,|2^. We 
have also investigated the distribution of delay times nu- 
merically and find the distributions of the Wigner delay 
time and the dwell time to coincide for energies far away 
from a band-edge for all disorder strengths. This, how- 
ever, breaks down for energies close to the band-edge and 
strong disorder, when the dispersive effects of the band 
structure deform the wavepacket so much so as to render 
the description in terms of the motion of a wavepacket 
meaningless. The concept of a dwell time, clocked by 
a counter such as the imaginary potential, however, re- 
mains meaningful even under such circumstances. 
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FIGURE CAPTIONS 

1. The distribution of delay and dwell times for reflection 
from a long disordered passive medium for wave energy 
at the middle of the band {E = 0.0, 1.0). 

2. The distribution of delay and dwell times for reflection 
from a long disordered passive medium for wave energy 
close to the edge of the band {E = 1.9, 1.99). 
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